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A novel Coulomb gas (CG) description of low energy QC D4,(Nc) is constructed. The construction 
is based on the dual transformation of the QCD effective Lagrangian. By considering a large gauge 
transformation, the charges of this statistical system are identified with magnetic monopoles which 
carry fractional charges of strength l/Nc- Furthermore, the creation operator which inserts the 
magnetic charge in the CG picture is explicitly constructed and demonstrated to have a non-zero 
vacuum expectation value, indicating that confinement does occur. The Wilson loop operator as well 
as the creation operator for the domain wall in CG representation is also constructed. Additional 
support for the D-brane picture suggested by Witten is also found. Lastly, the relation of the CG 
picture with the instanton-quarks is also discussed. 

PACS codes: 



I. INTRODUCTION 

Color confinement, spontaneous breaking of chiral symmetry, the U{1) problem, 6 dependence, and the classification 
of vacuum states are some of the most interesting questions in QCD. Unfortunately, the progress in our understanding 
of them is extremely slovif. At the end of the 1970s A. M. Polyakov j|] demonstrated color confinement in QED^; 
this was the first example in which nontrivial dynamics was a key ingredient in the solution. Soon after, 't Hooft 
and Mandelstam ^ suggested a qualitative picture of how confinement could occur in QCD4. The key point, the 't 
Hooft - Mandelstam approach, is the assumption that dynamical monopoles exist and Bose condense. Many papers 
have been written on this subject since the original formulation however, the main questions, such as, "What 
are these monopoles?" ; "How do they appear in the gauge theories without Higgs fields?" ; "How do they interact?" , 
were still not understood (for a recent review see |^]). Almost 20 years passed before the next important piece 
of the puzzle was solved Seiberg and Witten demonstrated that confinement occurs in SUSY QCD4 due to the 
condensation of monopoles much along the lines suggested many years ago by 't Hooft and Mandelstam. Furthermore, 
condensation of dyons together with oblique confinement for nonzero vacuum angle, 9, was also discovered in SUSY 
models (a phenomenon which was also argued to take place in ordinary QCD; see |^). In addition to forming concrete 
realizations of earlier ideas, the recent progress in SUSY models has introduced many new phenomena, such as the 
existence of domain walls |^ which connect two distinct 0- vacua. New insights into confinement was also recently 
given by Witten in which he argued that domain walls connecting two vacua labeled by k and k + 1 behave 
similarly to _D-branes on which the SQCD strings can end. It is tantalizing to suggest that such phenomenon also 
take place in QCD; indeed, in this paper it is argued that this in fact does occur. 

With such a wealth of information now available for SUSY gauge theories, it is interesting to ask how one can apply 
that knowledge to QCD. In a recent review Shifman ||] addressed certain aspects of this question and we refer to that 
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paper for further references and recent development. Thus far, most works have approached such questions by starting 
with a supersymmetric theory where some particular dynamical aspect of interest is already known; then introducing 
explicit supersymmetry breaking terms by giving the gluinos or squarks mass; finally, in the limit of infinite mass 
one recovers a non-supersymmetric theory with decoupled gluinos and squarks. Unfortunately, any calculations in 
this limit are beyond reach. Consequently, theories with a small SUSY breaking mass term are usually studied and a 
qualitative picture of what happens in QCD (and other non-supersymmetric theories) is extracted. This procedure has 
lead to success in the identification of the complicated vacuum structure for such theories: there are N inequivalent 
vacua for each given 9 ^ . These vacua are distinguished by the phase of the chiral condensate and there are domain 
walls which interpolate between them. It is believed that in QCD such a similar structure exists. 

Rather than beginning with the SUSY model and then breaking the symmetry, our analysis begins with non- 
supersymmetric QCD. For the purposes of the paper it will be assumed that the theory is in the confining phase and 
that chiral symmetry is broken. The question to be addressed is then: "Does the picture which emerges in SUSY 
models reflect what happens in QCD?" . We will demonstrate that the answer is indeed yes. This strategy clearly 
does not prove confinement, nor does it elaborate on the microscopic mechanism responsible for it, however, it will, 
however, demonstrate that the assumption of confinement implies the condensation of massless monopoles and vice 
versa. Furthermore, the role that the 9 parameter plays in the condensation of dyons with non-zero electric charge 
will be illuminated. In addition, this strategy will allow us to demonstrate that the states in QCD for each given 
0-angle is indeed classified by the integer number k labeling the phase of the dyon condensate. Finally, domain walls 
which interpolate between these different phases will be shown to have wall tension ~ N, and it is consistent with 
their interpretation as D-branes rather than as QCD string solitons. In addition, a relationship between the dual 
representation of low energy QCD, in terms of the Coulomb Gas of monopoles, and instantons will be discussed. 
This leads us to conjecture that at large distances our particles with fractional magnetic charges, from the CG 
representation, are in fact the instanton-quarks suspected long ago Q. 

Before going into a detailed analysis of the effect Lagrangian the starting point of our analysis will now be elaborated 
on. Experience in SUSY models demonstrates that the effective Lagrangian approach and duality transformations are 
very effective tools in the analysis of the large distance dynamics in the strong coupling regime (see e.g. reviews [^). 
These tools will be adopted for the study of QCD. The first key element in the analysis is the effective Lagrangian 
approach. As is well known, there are two different definitions of an effective Lagrangian in quantum field theory. 
One of them is the Wilsonian effective Lagrangian describing the low energy dynamics of the lightest particles in the 
theory. In QCD, this is implemented by effective chiral Lagrangians for the pseudoscalar mesons. Another type of 
the effective Lagrangian (potential) is defined as the Legendre transform of the generating functional for connected 
Greens functions. This object is useful in addressing questions about the vacuum properties of the theory in terms of 
vacuum expectation values (VEV's) of composite operators, as they minimize the effective action. Such an approach 
is well suited for studying the dependence of the vacuum state on external parameters, such as the light quark masses 
or the vacuum angle 9. However, it is not useful for studying, for instance, 5-matrix elements because the kinetic 
term cannot be recovered in such an approach. The utility of such an approach to gauge theories has been recognized 
long ago for supersymmetric models, where the anomalous effective potential has been found for both the pure gauge 
case 1^ and SQCD ||l^. Properties of the vacuum structure in SUSY models were correctly understood only after 
analyzing this kind of effective potential. The second key point of our approach is the observation that such an 
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effective potential in QCD (which is a collection of self-interacting Sine-Gordon (SG) fields) has a dual representation 
in terms of a Coulomb Gas of several species of charges, which will be described in detail in Section 3. 

Since our description starts from color-singlet fields (the phases of the chiral condensate), one might think that 
no information about confinement, which is clearly a color phenomena, can be extracted from such an analysis, 
fiowever, contrary to ones instinct, it is possible to obtain color information due to the existence of the free parameter 
9, which plays the role of the messenger between colorless and colorful objects. Indeed, the dependence of the 
physical observables as a function of 6 can be studied exclusively in terms of the colorless degrees of freedom by 
using the effective Lagrangian approach. Furthermore, the manner in which the theory transforms under both large 
gauge transformations and under the variation of the 9 parameter, when monopoles (colorful object) acquire a color 
electric charge, are known independently. Consequently, understanding the 9 dependence gives us an understanding 
of confinement, as 6 acts as a link between colorless and colorful degrees of freedom. 

This paper is an extended version of the letter and is organized as follows: In Section 2 we overview the properties 
of the effective Lagrangian of QCD Section 3 is devoted to derivation of the Coulomb gas representation of this 
effective Lagrangian, this representation can be interpreted as a dual form of the original low-energy effective field 
theory; the massless particles of the new statistical mechanical problem are identified with monopoles; the measure for 
this statistical ensemble is further analyzed and it is then argued that our particles can be identified with instanton- 
quarks. In Section 4 several consequences of the Coulomb gas representation are derived, in particular, the expression 
for the monopole creation-operator (magnetization) is obtained and its VEV calculated; The magnetic potential in 
the presence of the domain walls is also discussed and it is argued that the QCD string can end on the domain wall. 



II. EFFECTIVE LAGRANGIAN AND 6 DEPENDENCE IN QCD 

Our analysis begins with the effective low energy QCD action derived in |jl^, which allows the 6'-dependence of 
the ground state to be analyzed and is crucial to the identification of monopole charges. Within this approach, the 
Goldstone fields are described by the unitary matrix Uij, which correspond to the 75 phases of the chiral condensate: 
= -\^L-^n)\U,j with 

■7r"A" . 2 77' 



U = exp 



V2- 



In .JWf f^' 



UU+ = 1, (I) 



where A° are the Gell-Mann matrices of SU{Nf), tt" is the pseudoscalar octet, and = 133 MeV. In terms of U 
the low-energy effective potential is given by fl^ : 

Wqcd{0,U) = - Mm — logV exp IvE cos ( --9 + i-logDct U +— l] +-VTt (mU + m+U+)\ . (2) 
v^oo V [ \ P P p J 2 J 

All dimensional parameters in this potential are expressed in terms of the QCD vacuum condensates, and are well 
known numerically: m = diag(m* K^'^P') |); and the constant E is related to the QCD gluon condensate E — 
{bas / {32Tr)G'^) . The only unknown parameters in this construction are the integers p and q which play the same 
role as the discrete integer numbers classifying the vacuum states in SUSY theories. These numbers are related to 
a discrete symmetry which is a remnant of the anomaly, and can be found only by explicit dynamical calculations. 
Various arguments can be put forward which support different values of p and q, see e.g. [llSUsI; however, in what 
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follows p and q will not be fixed, rather only the ratio will be constrained so that in large Nc limit q/p ^ 1/Nc such 
that the U{1) problem is resolved. 

It is possible to argue that equation (^) represents the anomalous effective Lagrangian realizing broken conformal 
and chiral symmetries of QCD. The arguments are that: 
Equation (||): 

(i) correctly reproduces the WW effective chiral Lagrangian ||l^ in the large Nc limit 

[ For small values of (0 — ilogDet U) < Tr/q, the term with I — dominates the infinite volume limit. Expanding 
the cosine (this corresponds to the expansion in q/p ^ 1/Nc), we recover exactly the WW effective potential 
[0 together with the constant term —E = — (6as/(327r)G^) required by the conformal anomaly: 

Wvvw{0,U,U+) ^-E- ^^^^^{e~i\ogBetUf -^Tr{mU + m+U+) + ... . (3) 

here we used the fact that at large Nc, E{q/p)^ — —(y'^)YM is the topological susceptibility in pure YM theory. 
Corrections in XjNc stemming from Eq.(|^) constitute a new result of Ref. [p^.] 

(ii) reproduces the anomalous conformal and chiral Ward identities of QCD 

[ Let us check that the anomalous WFs in QCD are reproduced from Eq.(|2|). The anomalous chiral WI's are 
automatically satisfied with the substitution Q ^ {9 — i logDet U) for any Nc, in accord with |jl^. Furthermore, 
it can be seen that the anomalous conformal WFs of for zero momentum correlation functions of the operator 
in the chiral limit iriq are also satisfied when E is chosen as above. As another important example of 
WFs, the topological susceptibility in QCD near the chiral limit will be calculated from Eq.(^. For simplicity, 
the limit of SU{Nf) isospin symmetry with Nf light quarks, ^ ^qcd will be considered. For the vacuum 
energy for small < ir/q one obtains 

/ a \ ^ 



E.aci9) = -E + mg{^^)Nf cos i^— j + 0{m^) . (4) 

Differentiating this expression twice with respect to 9 reproduces the result of : 

Ihn ^ J dxe^^^O\T{^GG{x) ^GG(O)} |0) = -^^^ = ]^-.(**) + ^K) . (5) 

Other known anomalous WFs of QCD can be reproduced from Eq.(^ in a similar fashion. Consequently, Eq.(|^) 
reproduces the anomalous conformal and chiral Ward identities of QCD and gives the correct 9 dependence for 
small values of 9, and in this sense passes the test for it to be the effective anomalous potential for QCD. ] 

(iii) reproduces the known results for the dependence at small 9 , but may lead to a different behavior for large 
values 9>n/qifq^l 

[ As mentioned earlier, at small 9 ^ ir/q our results are identical to those found in [Q; the main difference 
with ||lj] arises when 9 — ilogDetC/ ^ 1 where cusp singularities occur. These singularities are analogous to 
the ones arising in SUSY models and show the non-analyticity of the 9 dependence at certain values of 9. The 
origin of this non-analyticity is clear, it appears when the topological charge quantization is imposed explicitly 
at the effective Lagrangian level. Thus, the cusp structure of the effective potential seems to be an unavoidable 
consequence of the topological charge quantization (which was not explicitly imposed in the approach of ] 
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An interesting note is that in general the 9 dependence appears in the combination 9 /Nc which naively does not 
provide the desired 27r periodicity for the physical observables; however, such a behaviour is derived from Eq.(|2|) which 
does have 27r periodicity. This seeming contradiction is resolved by noting that in the thermodynamic limit, ^ oo, 
only the term of lowest energy in the summation over I is retained for a particular value of 0, creating the illusion 
of 9/Nc periodicity in observables. Of course, the values 9 and 9 + 2tt are physically equivalent for the entire set of 
states, but not for a selected individual vacuum state. Consequently, the 9 dependence in the infinite volume limit 
appearing in the combination 9 /Nc is a result of being stuck in a particular state. The reader is referred to the original 
papers for more detailed discussions of the properties of the effective potential (||). One final point to be made 
before moving on is that in general Eyac{9) is a multi- valued function with cusp singularities aX 9 — (^''+^)^ where 
one solution switches to another. Furthermore, the number of additional meta-stable vacuum states is very sensitive 
to the integer parameters p and q introduced earlier. In particular, for the physical values of the quark masses, there 
are q — 1 additional local minima of the effective chiral potential, which are separated by large barriers of strength 

E from the true physical vacuum of lowest energy, but which are almost degenerate in energy, Ai? m. Such 



availability of almost degenerate vacua imply the appearance of QCD domain walls and was studied in |17|. The 
appearance of these walls have interesting consequences in this work as well. 

This section is closed with a few technical but important remarks: The effective low energy Lagrangian described 
above has a very special SG structure. The fact that such a structure appears in the second term of Eq. (H), 
m{U +U^) ^ "^i cos(0i), is quite natural and is associated with the Goldstone origin of the (pi fields. The first term 
has a similar SG structure ~ cos 4>i ^ ^) /^c)', however, its origin is less trivial and requires some explanation. 
Firstly, WI's imply that the singlet combination J^i'^'i always appears with 9 in the form {9 — ilogDctJJ) [Q. 
Secondly, the appearance of the cosine interaction, cos{9/Nc), leads to the following scenario in pure gluodynamics 
((^i's are frozen): the (2fc)'^ derivative of the vacuum energy with respect to 0, as — s- 0, is expressed solely in terms 
of one parameter, for arbitrary k: 



QQ2k 



?=0 



J l[d'x, {Q{x,)...Q{x2k)) - [^J , (6) 



where, Q ^ Gf^iyG^^. This property was seen as a consequence of Veneziano's solution of the U{1) problem [[l8[ . 
The reason that only one factor appears in Veneziano's calculation is that the corresponding correlation function, 
^ J Y\^!Lid''^Xi{Q{xi)...Q{x2k)) , becomes saturated at large distances by the Veneziano ghosts whose contributions 
factorize exactly, and was subsequently interpreted as a manifestation of the 9/Nc dependence in gluodynamics at 
small 9. However, at that time it was incorrectly assumed that such a dependence indicates that the periodicity in 
9 is proportional to Nc- Later on it was argued that the behavior in Eq.(^) is a consequence of the holomorphic 
structure for the non-perturbative part of the QCD partition function where the non-perturbative vacuum energy 
depends only on a single complex combination r = I/^q + «3^y in the same way as it was in SUSY modelsl^ The 



'^There is an essential difference between SUSY models and QCD however: in SUSY models holomorphy is an exact property 
of the effective super-potential; however, in QCD it is only a property of the non-perturbative effective large distance part, see 
H, for more details. 
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behaviour in (|6|) was also demonstrated by "integrating in" a very heavy fermion and integrating it out afterwards, 
which must leave the gluodynamics unaltered All three independent arguments support the property (|^), 

consequently Eq. (Q) will be used as the defining effective action for low energy QCD. 



III. COULOMB GAS REPRESENTATION OF LOW ENERGY DYNAMICS 

The effective low energy QCD action ^ has a very special quality, all the interactions have a trigonometric form 
and are in particular cosinusodal. This is the defining character of a SG model. The presence of many fields and 
many cosine terms with different harmonics which interact in a highly non-trivial manner serve to slightly complicate 
the situation. Nevertheless, many of the special properties of the SG theory apply to this model, the admittance of 
a Coulomb gas (CG) representation for the partition function is no different. Although this is a four dimensional 
theory, and questions about renormalizability of the theory may come to mind, there are no such issues here since 
the effective action is a low energy one. Following the usual procedure for mapping a statistical CG model into the 
field theoretic SG model, the CG picture that arises from the effective low energy QCD action, Eq. (||), will be 
derived in this section. The statistical model will be seen to contain several species of charges which appear due to 
the presence of several cosine interactions in the field theory model. The physical meaning of these charges will be 
illuminated by applying large gauge transformations on the partition function and determining the manner in which 
the charges transform. This will lead to the identification of one of the charge species as magnetic monopoles and 
of the singlet combination, 0, in the SG model as the magnetic scalar potential. The remaining charges will also be 
argued to have magnetic properties; however, it will prove difficult to construct a precise statement about them. The 
statistical model will then be analyzed in subsequent sections of this paper. 



A. Formal Derivation 



Although the mapping between a SG theory and its CG representation is well known 0, in this section several 
illustrative steps in its derivation will be given for completeness. The existence of many fields and cosine terms only 
serve to make the formulae more bulky; however, the basic strategy is the same as in the standard case. Using the 
effective potential in Eq. (j^) the defining partition function is taken to be, 

Z — cxp I ^ y d'^x^Ecos. {(j> — 9 + 2TTn) + 27rfc^ + mi cos{4>i) + . . . + ruN; cos((/)Ar^)|| ^ (7) 

Here, the matrix U characterizing the phases of the chiral condensate has been placed in diagonal formQ, U = 
diag((/)i, . . . , (f>Nf ), the combination of fields (p ~ Tr U represents the singlet, — m°|(5'°^'°)| denotes the mass of 



^It is well known that this is the most general form for i7-matrix describing the ground state of the system jlj]. The off- 
diagonal elements of U describe the fluctuations of the physical Goldstones which are neglected. This is allowed since only the 
diagonal elements are relevant in the description of the ground state which is the focus of this work. This truncation can be 
justified a posteriori by demonstrating that the classification of vacuum states based on the CG representation exactly coincides 
with the classification based on the effective Lagrangian approach, see Eq. where only the diagonal elements are relevant. 
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the a-th flavor of quark together with its condensate, E = (6as/(327r)G^) is the vacuum energy of the system, and 
the angled brackets represent the Gaussian weighting, 

(■••)=E E h'^'^ ■■■ e-^' I (8) 

fc=0 n= — oo 

where is a coupling constant ~ so that the fields (f>a are dimensionless and all distances are measured in the 
units of f^^. 

The following strategy for obtaining the CG representation from the SG model will be employed here: a series 
expansion in the vacuum energy and the quark masses will be carried; the cosine terms will then appear outside of 
the exponential and are weighted by the Gaussian weighting; introducing ^2 valued fields to represent the cosine 
terms allows the dynamical fields to re-appear in the exponential, however, now they appear linearly; in this form it 
is possible to completely integrate out the dynamical fields which leads to the final CG representation. These steps 
are now carried out in order. Firstly, performing a series expansion in E and rua independently leads to the following 
form, 

1 

{Mo,...,MjVj.}=0 



Y[(^ldx<-^^cos{U^^^^))Y^\ (9) 



a=l 

By introducing the 2^2 valued fields Q^°'\ for a — 0, . . . ,Nf, to replace the cosine interactions the dynamical fields 
can be integrated out exactly. For a 7^ the fields are introduced via. 



Q(»)=±l 

and the charge species (3*^°^°^ will be said to be dual to the field cj^a- The singlet case, a = 0, contains additional 
parameters, 

cos (^4, - 6 + 2^n) + 2^k \ = ^ E e^'^'^'^^-^+^-^^+^^l'^). (10) 

and the charge species Q'^'^^ will be said to be dual to the singlet combination, (j) — Tr U . The explicit presence of the 
0-angle in ( [l0| ) is essential in identifying the physical meaning of this charge species. Inserting both expressions into 
the series expansion of the partition function, Eq. (^), leads to a form in which the dynamical fields appear at most 
quadratic in the exponential and they can be completely integrated out. 



Furthermore, as will be demonstrated shortly, the most important contribution is related to the singlet field (j) — Tr U which 
is unambiguously defined. All contributions related to the non-singlet fields are suppressed in the chiral limit by the quark 
masses rrii. 
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N 



f 



{Qf ,..'.q1"^'} = ±1 

The functional integral is trivial to perform and one arrives at the dual CG action. The brackets can then be removed, 
as long as the summations over n and k are restored. After some rewriting the CG action is, 

\ n=-oo / ^ a=l I b,c=l 6=1 c=l 

+ ^Qi'^)G(4")-4"))g('^) . (12) 

fc,c=l J 

Here, Q^"* = X^b^i QI""* ^he total Q^^^ charge for that configuration, G{x — y) denotes the relevant Greens function 
of the Laplace operator (— □G(x — y) — 5{x — y)), and the factor of p appears due to the redundancy introduced by 
the summation over k in Eq.(^). 

Notice that the fugacities of the charge species Q^"^"^ are given by the masses of the a*** quark, while the fugacity 
of the charge species Q^"^ is proportional to the gluon condensate E. An important point is that the fugacity of charge 
species Q^"^^^ vanishes in the chiral limit, while that of the Q'"^ charges remains non-zero. A second point is that, in 
the chiral limit, this representation does not obviously have invariance in 9, while in QCD the 0-angle appears with 
quark masses and hence disappears in this limit. It is possible to demonstrate that the partition function has this 
property; however, because its proof is rather technical, the details are deferred to appendix A. 

There are several important features of the action (112) which should be noted. Firstly, the summation over n forces 



the total Qf^"^ charge, to be an integer. Such a constraint is the analog of the quantization of the topological 

charge and is to be expected. Additionally, since this species has charges ^ p ^ V^c, this constraint enforces a 
fractional quantization on the total charges: the difference in the number of positive and negative charges of species 
(5^°^ must be an integer multiple of Nc- Secondly, the 6'-angle only interacts with the Q(°^ species. This directly 
results from the species being dual to the singlet field 0, which is the only field directly interacting with the 
6'-angle. Finally, the 6'-dependence acts only to supply an overall phase factor for each configuration and leads to the 
very natural interpretation of non-trivial ^-angles as introducing an overall background charge. Turning attention 
to the interactions amongst the charges, the species Q'-^-' is seen to interact with all species Q^"'' (« = 0, . . . ,iV/); 
however, the other species, Q^°'^°\ only interact with their own species and species Q'^^\ This peculiar behave is, once 
again, due to species Q^°^'s association with the singlet while the other charge species are associated with particular 
components of the chiral condensate. It is interesting that the structure of interactions between the charges resembles 
the 't Hooft "determinant" interaction, in which all Nf fermions form a vortex with Nf legs, while different species 
do not notice one another. Such a structure can be easily understood in the instanton picture, in which the existence 
of exactly one zero mode for each given flavor provides precisely this kind of interaction. In the present context, 
however, a straightforward explanation in terms of the original gauge and matter degrees of freedom is not available. 
Nevertheless, the interactions occurring in ( |l2| ) are quite similar to other known statistical systems, yet they do not 
literally correspond to interactions appearing in any previously discussed model. 

Figure |l| summarizes the manner in which interactions take place in our statistical system. Notice that the SU{Nf) 
symmetry is restored in the limit of equal quark masses. This occurs since the fugacities of all charge species are 
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identical in this limit, and the multiplicity of any of the diagrams depends only on the number of legs and not on 
what charges are attached to them. 




FIG. 1. The interactions amongst the charges in the CG representation. In this diagram a 7^ fo 7^ c 7^ 0. The charges can 
interact only amongst themselves (diagrams furthest left), or with the Q'"-* charge. 



B. Physical Interpretation of Charges 

Expression (^2|) clearly shows that the statistical ensemble of particles interact according to the Coulomb law 
\x — y\^^ in static (time-independent) configurations { ^ \x — y\^'^ for time varying ones). An immediate suspicion 
following from this observation is that these particles carry a magnetic and/or electric charge, since charges of that 
type interact precisely in the above manner. Due to the direct interaction of the 6'-angle and the charges, Q^"-*, 
associated with the singlet field, (f>, it is quite plausible that they may in fact have a magnetic nature. This suspicion 
will be corroborated in a moment. Furthermore, since there are several species of charges in the CG picture, a 
natural question is, "What is the physical relevance of the remaining flavor non-singlet charges Q^"'^"'?" Evidence of 
their magnetic nature is found by noting that the VEV's, {(paiO)), of these fields do in fact depend on the 0-angle. 
The implication is that all charge species, not only the one dual to the singlet field, have a magnetic origin and are 
connected to monopoles. Unfortunately, a precise relation of the charges dual to the non-singlet fields and monopole 
charges is still missing. 

As alluded to above, a precise relationship of the Q'^^^ charges and magnetic monopoles can be made. However, in 
order to prevent a long interruption in the presentation, only the simplest identifications will be given here, while the 
interested reader is referred to Appendix B for further detailed discussions. The charge was originally introduced 
in a very formal manner so that the QCD effective low energy Lagrangian (|^) can be written in the dual CG form 
(|l^). Now the physical content behind these manipulations will be given. 

To begin, consider the SU{2) Georgi-Glashow model in the weak coupling regime, with a 6'-term, and let the scalar 
(f>° have a large VEV. The monopole solution can be constructed explicitly and the so-called Witten effect , where 
the monopole acquires an electrical charge, takes place^ Let N denote the generator of large gauge transformations 



^It should be noted that the following arguments hold not only for the Georgi-Glashow model. Monopoles could appear in 
the system via a nontrivial dynamical effect and are not necessarily described in the terms of the original fields in the weak 
coupling regime. For example, they can appear as a result of gauge fixing which is certainly not the unique gauge-invariant 
construction. Nevertheless, it is believed that the Witten effect takes place for such 't Hooft monopoles as well as for any 
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corresponding to rotations in the U{1) subgroup of SU{2) picked out by the gauge field, i.e. rotations in SU{2) about 
the axis n'^ = Rotations by an angle of 27r about this axis must yield the identity for arbitrary configurations, 

which implies ||l9| that the magnetic monopoles carry an electric charge proportional to 6. Indeed, 

l = e'2.jv^g.2.f-.e^^ (13) 

where, 

are the magnetic and electric charge operators respectively, expressed in terms of the original fields, and v is the VEV 
of ^'^ at infinity. The combination ^ — n„j in Eq. (|l^) is an integer and determines the magnetic charge of the 
configuration. As usual, it is assumed that ( |l3|) remains correct in the strong coupling regime when v is not large 
and/or in the more radical case when <i>° is not present in the original formulation. Indeed, as explained in |^ the 
existence of ^'^ is not essential and some effective fields may play its role. One finds that monopoles do exist and the 
Witten effect expressed by formula ( p^ ) remains unaltered even when monopoles appear as singularities in the course 
of the gauge fixing procedure as described in 

Restricting attention to terms which are proportional to the ^-parameter, a comparison between the CG represen- 
tation, Eq.(^2|), and Eq.(|l3|) will now be carried out. From the CG the relevant term is the total charge, of the 
configuration, while in Eq.(^3|) the relevant factor is the total magnetic charge ^ for each time slice. The following 
identification is then made, 

Qt^^= "™ € (15) 
47r 

A non-trivial check on this identification can be made by noticing that the Q^^^ species are fractionally charged 
q(o) ^ 1/Nc, while the total charge Q^-'is forced to be integer. This follows from the constraint that the difference 
in the number of positive and negative charges of species Q must be an integer multiple of p ~ Nc (see discussions 
after (|l^)) and this behaviour is identical to the that of the total magnetic charge. 

From these simple observations one can immediately deduce that our fractional magnetic charges Q^^'' cannot be 
related to any semi-classical solutions, which can carry only integer charges; rather, configurations with fractional 
magnetic charges should have pure quantum origin. Of course, this is a simplified explanation of the identification 
between the charges Q'^*'-' from ( [l^ ) and the physical magnetic charges. To make this correspondence more precise the 
transformation properties of each relevant degree of freedom under a large gauge transformation must be computed 
and the details of such a calculation can be found in Appendix B. 

As already mentioned, the charges Q^"'^"-' are also suspected to carry magnetic charge, however, a precise statement 
cannot be made because the large gauge transformations are sensitive only to the singlet combination of fields and 
consequently only to the dual charge (5*-°^. 



other monopoles. The most important feature of the Witten effect are the properties of large gauge transformations rather 
than specific features of the theory and matter contents. 
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It is interesting to note that a similar phenomenon in the two-dimensional 0(3) cr- model (more generally, in CP^ 
models) has been known for a long time |20[ . Namely, an exact accounting and resummation of the n-instanton 
solutions maps the original problem to a 2d-CG with fractional charges (the so-called instanton-quarks) ; however, the 
total topological charge of each configuration is always integer. In that case, it is clear that the quantum fluctuations 
completely reconstruct all degrees of freedom: each instanton is essentially a superposition of Nc instanton-quarks 
such that, locally, they appear as fractionally charged particles. This superposition is quite nontrivial; however, it is 
stable under small fluctuations and eventually becomes the only relevant degree of freedom. Indeed, the corresponding 
two-dimensional CG of the instanton-quarks can be rewritten as a quantum field theory with effective Lagrangian 
given by the SG model where the dynamics is carried by some effective field <j) |20j. This effective SG Lagrangian 
describes confinement (mass-gap) and many other important properties of the 0(3) cr-model quite well. Our CG 
system with fractional magnetic charges, Q'^^-' ^ 1/p ^ is very similar to the situation described above. The 

difference with is the starting point. In the case of the 0(3) a-model the statistical-mechanical problem was 
reduced to some low-energy Lagrangian. However, in the present case we start from the low-energy Lagrangian and 
a statistical-mechanical model of fractional magnetic monopoles is derived. 

Until now, the integers p and q have been considered as free parameters in the theory, and can only be fixed by 
an explicit dynamical calculation. The only constraint on these parameters has been that in the large Nc limit their 
ratio satisfies q/p ^ 1/Nc so that the U{1) problem is resolved. An argument which supports a particular choice for 
these arguments will now be given. Recall that q/p is related to the 9 dependence of physical observables (such as the 
vacuum energy Eyac{d) as can be seen from Eq.(^)). In particular, an observable with expectation value ~ exp(i-j|-) 
corresponds to q = 1, p = Nc. The analysis of softly broken SUSY theories |^ strongly suggests that this choice of 
parameters, q = 1 and p = Nc, is in fact correct even for non-supersymmetric models. In this case the number of 
integrations over d'^xf^'^ in Eq. (pT|) exactly equals ANck, where k is an arbitrary integer (not to be confused with the 
summation parameter appearing in the partition function) . This is a direct consequence of the fact that the number 
of species Q'*'-' must be an integer multiple of N^ This number, ANck exactly corresponds to the number of zero 
modes in the fc- instanton background l2l| ], and, correspondingly, to a number of collective variables associated with 



*To be more precise, the difference in the number of positive and negative charges of species Q^"' must be an integer multiple 
of Nc- However, in what follows charges with a specific sign, say positive, will be identified with instanton-quarks. Such an 
identification will be made by only retaining a single term in the expansion of the cosine function in (^o|) to obtain a statistical 
model with only like signed charges analogous to (jll| [l2| ). In this case the number of particles, and not the difference in 
the number of positive and negative charges, must be an integer multiple of Nc- Restricting the model to contain only like 
charged particles is slightly illegal for = in the presence of light quarks, because in the large volume limit the corresponding 
contribution to the partition function is suppressed by the inverse volume of the system (due to the neutrality requirement) 
and/or by the quark mass. However, the question regarding a measure of the corresponding configuration (before integration 
over d^xf^) is a perfectly acceptable one which has an answer. In general, it is expected that both instantons and anti-instantons 
must be present in the system in order to give a nonzero result for the partition function in the chiral limit (see discussions 
at the end of this section). In this respect there is a difference with the 0(3) a-model without fermions where exclusively 
instantons could provide a nonzero contribution to the partition function. 
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these zero modes. In other words, 4Nck coordinates x\ describe the collective coordinate integration measure of the 
fc-instanton solution. 

Motivated by the analysis of the two-dimensional ct- model in |2^] and by the above observation regarding the number 
of collective variables AN^k mentioned above, we conjecture that at large distances our particles, Q^^\ with fractional 
magnetic (and electric for 9 ^ Q) charges are indeed the instanton-quarks (i.e. these are related to instantons) 
suspected long ago Q]. One immediate objection to this conjecture is that since it has long been known (see e.g. 
p2[) that instantons can explain most low energy QCD phenomenology (chiral symmetry breaking, resolution of the 
U{1) problem, spectrum, etc) with the exception of its most important property -confinement; and we claim that 
our magnetic monopoles, Q^"-*, are instanton-quarks and yet also claim that confinement arises in this picture (they 
will Bose-condense, see section IV); how can this be consistent? This seeming objection is answered by noting that 
it in the dilute gas approximation, when the instantons and anti-instantons are well separated and maintain their 
individual properties (sizes, positions, orientations), quark confinement can not be described. However, the lessons 
from the two-dimensional a model teaches us that in strongly coupled theories the instantons and anti-instantons 
lose their individual properties (instantons will "melt" ) their sizes become very large and they overlap. If this happens, 
the description in terms of the instantons and anti-instantons is not appropriate any more, and alternative degrees 
of freedom should be used to describe the physics. The relevant description is that of instanton-quarks (which, 
according to our conjecture, are particles with monopole charges Q*-*^-*). Further to this point, lattice simulations do 
not contradict this picture where large instantons induce the magnetic monopole loops forming large clusters, see 
e.g. and references therein. Also, the connection between monopoles and instantons on the classical level is not 
a very new idea ||2^ . Indeed, quite recently, such a relation was established for the periodic instantons (also called 
calorons) defined on x in the presence of a Wilson line p5|| . Furthermore, a similar relation was seen in the 
study of Abelian projection for instantons ||2^,^, albeit at the classical level. In particular, Brower et. al. [2^ ] 
demonstrated that the instanton's topological charge, Q, is given in terms of the charge, M, of the monopole, which 
forms the monopole loop, by the expression Q = This formula is very similar to our relation (p^), where the total 
topological charge, Q^j^\ for a configuration containing a number of particles, {Q^^"*}, described by the system ( [T^ ) 
was identified with the total magnetic charge for each time slice for the same configuration. In spite of this similarity, 
the physical interpretation of this relation is quite different for these two cases. In the former case it is an identity 
for a configuration with integer monopole and topological charges satisfying the classical equations of motion; while 
in the latter case it is a relation for a configuration of the fractionally charged monopoles and topological charges 
(fractional) for the instanton quarks (which presumably appear after integration over quantum fluctuations in the 
multi-instanton background). However, given that the instantons at the classical level, and the instanton-quarks at 
the quantum level both have similar relations suggests that our conjecture is possibly correct. Additional supporting 
arguments based on the instanton measure for arbitrary gauge groups are given in our concluding remarks. 

As a last remark regarding our conjecture, note that the QCD effective Lagrangian (^, at = 0, is invariant 
under CP transformations; consequently, the CG representation ( pi] ) is also invariant under CP a± — 0. This 
implies that only when instantons and anti-instantons are taken into account can the statistical ensemble of charges 
be equivalent, at large distances, to a statistical ensemble of instanton-quarks. If either are left out of the description 
CP symmetry is broken. The same conclusion also follows from the fact that the partition function, (|l^), is non-zero 
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in the chiral Umit^ nia = . Such a behaviour of the partition function imphes that our system (yjj) can be equivalent 
to the statistical ensemble of the instanton-quarks if and only if both instantons and anti-instantons are present in 
the system, and there is no excess of the topological charge (for 9 = 0) which, otherwise, would lead to the chiral 
suppression of the partition hmction. 



IV. OPERATORS 



In the previous section, a CG representation which describes the low energy dynamics of the QCD action was 
derived. The charges were found to have magnetic properties, and the charge species, Q'-*'-', dual to the singlet field, 
(j), was explicitly identified as magnetic monopoles. However, the formal mapping from the SG description to the CG 
description, to (|l^), is not very useful if the dynamics of the fields and charges are not understood. In order to 
understand the long-distance features of the theory, the manner in which the various charge species behave must be 
investigated. In particular, if the magnetic charges Bose-condense, this indicates the onset of quark confinement. To 
investigate the possibility for such a condensation an expression for the magnetic charge creation operator, M, must 
be found and its VEV (magnetization) calculated. In this section the magnetization (Ai) will be demonstrated to 
be non-vanishing; consequently, the theory is in the confining phase. Of course, this is not a proof of confinement in 
QCD because the effective action was constructed under the assumption that the system is in the confining phase, 
i.e. (^) contains only colorless degrees of freedom. In spite of this, the result is quite nontrivial and can be considered 
as a self-consistency check of our identifications. In addition, it provides a very nice intuitive picture of the origin of 
the confinement in QCD. 



A. Charge Creation 

Consider inserting a charge, qo = of species Q^"^ in the bulk at the point X. To accommodate this insertion, 
the partition function must be altered by restricting, say, the coordinate Xq"'' = X and q'^^ = qo. Tracing the steps in 
the previous section in reverse order with this restriction leads to the appropriate operator equivalence. However, it 
is simpler state the result and then demonstrate that it is correct. The operator equivalence, which is proven correct 
shortly, is, 

M(qo,X) =e"i"(*(^)-''+'""+2"i'=), (16) 

where M denotes the operator which inserts the charge in the bulk. The singlet combination naturally appears here 
since (5^°^ is dual to (f>, while the appearance of 9, n and k appear due to the manner in which Q^o) was introduced 
originally (see Eq. (|lo|)). Clearly this operator must be inserted under the summations over n and k. On inserting 
( p^ ) into the partition function and performing the series expansion in E and ma then introducing the 7L2 valued 
fields, as before, the modified CG action is arrived at. 



^However, only one species of particles with charges Q*^"' contribute to the partition function in tills case; other species with 
charges Q^"^"-* have zero fugacities in the chiral limit and, therefore, are not present in the system. 
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The term appearing outside the inner parenthesis is the suggested operator for the creation of monopoles. Since X 
and qo are not being summed over, this term wiU act as a source for a fixed charge at a fixed position. Collecting 
this term with the first term and performing the functional integrals over the scalar fields, demonstrates that qo will 
interact within the statistical model precisely in the form of charge species Q*^"-* with the exception that its charge 
and position is not being summed over. Consequently, the operator relation ( p^ ) is proven to hold. 

The insertion of any other charge species Q^'^^^^ can be found via an analogous method, and the operators for these 
charge species are found to be, 

X(q„X) ^e'l--^-^^) (17) 

Notice that ( [T^ ) explicitly depends on 9, n, k and the singlet combination (/>, while ( |T7| ) only depends on the scalar 
field of its own species, 4>a^ and not on any other parameters in the theory. This is of course a consequence of Q^*'-' being 
associated with monopoles and must transform correspondingly under large gauge transformations. Furthermore, the 



manner in which the singlet field, 0, appears in (16) allows the further identification of 4> in the CG representation as 
the magnetic scalar potential. 

Now that the relevant operator relations have been identified, questions about the dynamical nature of the system 
can be asked. The most pressing one being, "Does the operator which creates a charge have a non-zero expectation 
value?" This is a difficult question to answer in the CG representation, however, in the dual SG representation an 
answer at the semi-classical level can be given. Its solution, however, requires a knowledge of the VEV's of the phases 
of the chiral condensates, (0a), such that at the semi-classical level {M.) ~ e*'i"^'^"\ and the problem is reduced to 
the calculation of {(t>a)- These VEV's can be found by the minimization of the effective potential (^, and are given 
by the equations |l^ , 

sin(^((0)-^?) =-^^sin(0,), 0^TrC/-V0a, a=\,...Nf (18) 

where the potential has been chosen to be in the particular branch fc = 0. In the limit of equal quark masses, = to, 
and assuming ^ ^ 1, the system can be solved perturbatively, and to lowest order the VEV's are: ((/)a) — for 
small 0. Earlier on it was mentioned that the VEV's of all components of the chiral condensate acquire dependence 
on the ^-angle and this is now explicitly demonstrated. The suspicion of the magnetic origin of the non-singlet fields 
can now be justified somewhat, although an explicit construction, as in the case of the singlet field, is still missing. By 
inserting these VEV's into the expressions for the charge creation operators the magnetization at the semi-classical 
level can be evaluated. 



(Al(qa,X))^ <^ , (19) 




From Eq. (|19|) one can clearly see that operator M. is the order parameter of the system and its nonzero VEV implies 
the condensation of monopoles (dyons for non-zero 9). Eventually this condensation leads to confinement (oblique 
confinement for non-zero 0) in the system. 
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To close this section note that the phase of the order parameter (A^) (19) which labels different vacua in the 
CG representation (|l^) is determined by the same equation ( |l8|) describing the chiral condensate phases in SG 
representation (Q), which is the conventional description of low-energy QCD in terms of the chiral condensates. The 
physical meaning of the phases, (/)o, presented in the CG representation ( pT| ) and in the effective Lagrangian approach 
is, however, quite different. In the former case they represent the magnetic scalar potential; while in the latter 
case they are the phases of the chiral condensate. It is quite remarkable that in spite of the very different physical 
interpretations suggested by the two (dual) approaches lead to the same classification of vacua. One final remark is 
that generalizations of our calculations of (A^(qa,X)) for the case of different branches, arbitrary 9 and non-equal 
quark masses is quite straightforward, and reduces to the solution of Eq. ( p^ for the phases (0a)- The corresponding 
analysis has been carried out in ||l^ and it is not repeated here. 

B. Background Fields and Wilson Loop Operator 

As demonstrated in the previous section, the VEV of magnetization is non-zero, (A^(qo,X)) ^ 0; therefore, our 
system is in the confinement phase as expected. Consequently, the VEV of the Wilson loop {W) must show an area 
law dependence, and it is interesting to find an explicit form of the corresponding Wilson loop insertion in the SG or 
CG representation. Naively this seems like a hopeless goal, because the starting point, the effective low-energy QCD 
Lagrangian, does not contain any fundamental degrees of freedom (gluons and quarks) in terms of which the Wilson 
operator is defined. Nevertheless, as will be seen in a moment, the corresponding Wilson loop insertion operator can 
be recovered. The key point in this correspondence is again the identification of our charges Q^^^ with the physical 



magnetic charges, Eq. (|T5|). Once this identification is made, the singlet field, 0, in the expression (16) for the 
magnetic creation operator is interpreted as the magnetic scalar potential such that ^^(qojX) = e'*"^^-^^ represents 
the interaction energy of a monopole of (fractional) charge go &t position X in the presence of the potential of other 
monopoles. With this identification it is quite obvious what kind of replacement in our formulae (|^) and (^ij) should 
be made in order to insert a Wilson loop operator. Physics suggests the following picture: For each given time-slice, 
xq, the Wilson loop inserted (at {xi, X2) plane) into the monopole plasma behaves like a sheet of magnetic dipoles, and 
produces at large distance from the sheet a dipole magnetic field. The monopoles of the vacuum plasma, however, are 
polarized by this dipole field and react to produce a dipole field to try to cancel the field from the sheet. Therefore, 
this interaction of between a source with magnetic charges from the system changes the magnetic scalar potential (p. 
This change, however occurs only in the thin region where the magnetic plasma of the system does not cancel the 
field from the sheet. The thickness of this region is order , the only dimensional relevant parameter in the system. 
Therefore, the Wilson loop insertion changes the magnetic scalar potential, </>, by -I- where r] is an addition to the 
magnetic scalar potential due to a dipole layer of unit strength on the sheet S represented by the Wilson surface. 
The key point is that even though gauge fields are not present in our description, through the identifications of the 
charges as magnetic monopoles, just enough information about them is encoded in the magnetic scalar potential, 0, to 
obtain interesting properties. The transformation property of cf) under the Wilson loop insertion, as explained above, 
la (j) ^ (jy + rj. The function ri{x) has a discontinuity of 47r when the point x crosses the surface bounded by the Wilson 
loop for each given time slice. In a real physical situation the dipole layer is not infinitely thin (as mentioned above 
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it is order of E~^); however, incorporating the finite thickness is beyond the present analysis. 

An exphcit representation for the appropriate source term will now be constructed. It will be explicitly demonstrated 
that such a source interacts with the magnetic charges (5*-°^ exactly in the same manner as the <j) field, ~ e*'^*"'''. 
Therefore, the rj source insertion can be identified with an appropriate change in the magnetic scalar potential 
(j) ^ (p + rj. But, as explained above, the Wilson loop insertion must also lead to precisely this kind of shift <j) ^ (p + ij 
in SG (0) or CG (|ll]) representation. 

The insertion of background dependence is furnished by altering the Gaussian weighting (|^) to include source terms, 

p— 1 oo 



E E 

fc=0 n= — oo 



The classical equations of motion are of course modified by this insertion, and with appropriate choices of the source 
fields, 'r]^\ the classical solution can be made kink like, for example. However, rj^a"^ will not be constrained at this 
moment, and any source (Wilson loop, domain wall or anything else) can be described by this modification. Of 
course, it is interesting to study what this modification implies for the CG representation. On expanding the partition 
function in powers of E and nia, and introducing the ^2 valued charge fields one arrives at precisely equation ( [Tl| ) 
where the angled brackets are the modified ones above. Performing the functional integral posses no difficulty and 
one arrives at the CG representation, 

Z- V (m^nr^ {m^,/2r^f f / (o) (o) A / (^,) 

{Mo,...,Af„^.}=0 U i J^f J 

Nf 



a=l 



{Q['\...d"f'}=±i 



In the above, Scg is the CG action containing the interactions of the charges and is unaltered from its form given in 
equation (p^, also rj = rji + . . . + t/Nj. This then allows the identification of the operator, in the CG representation, 
which introduces a source for the classical equations of motion of the phases of the chiral condensates, 

Nf 

Cirj^)^e~'^I''^ ^^^^^ e'Er^.Q^-'-'C^"') J|e^E^iQi"'''»(4"') (20) 

in this formula the term e~~ J ^ naOri^ j-gpi-gggj^^g ^ self- interaction of the source with itself, it is irrelevant for the 
present purposes and will be omitted in the analysis. The second term in ( |20| ) is much more interesting and gives 
a very simple prescription for the insertion of any source. If one wants to insert some classical field configuration 
(source) in the original formulation, then in the CG picture one need only introduce a phase factor which gives rise 
to a source term for that background - all interactions among the charges remain unaltered. The charges simply 
interact with the source field just as if it were an external magnetic field, each charges species, Q'-°'^'^\ interacts with 
the source component with which it is associated to, r/a, while the singlet charge species Q^"-* interacts with all of the 
background fields. The special form of the SG action is what admits this very simple picture and is at the heart of 
the CG representation itself. Consequently, it is expected that the analysis carried out here has wide applicability 
and can be applied for any kind of sources. 
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In deriving (20) no assumptions about the source fields rja have been made. The remainder of this section will 
be devoted to the most interesting source - the Wilson loop operator insertion. As argued above, any source can be 
accounted for by the shift (j) (j^ + Tj, and the Wilson loop operator insertion is no exception. Our problem, therefore, 
is reduced to the calculation of a difference in magnetic scalar potential due to the insertion of this specific source. 
Before dealing with the Ad case, some relevant formulae from the 3d case, where the physics is well understood, will 
be reviewed; then the appropriate generalizations to the 4c? case will be made. The reason that the 3d analysis can 
be easily generalized to the 4d case is that in both cases the interactions between particles in the CG representation 
( |l^ ) is determined by the relevant Greens function {l/\x — y\ or l/{x — yY respectively); consequently, the Wilson 
loop insertion in the CG representation should be expressible in terms of the appropriate Greens function. Hence, if 
the Wilson loop insertion operator in 3d can be expressed exclusively in terms of the 3c? Greens function, l/\x — ?/|, 
the transition from three to four dimensions will be furnished by the replacement of the 3d Greens function by the 
4d Greens function such that the static limit is reproduced. In addition, it is expected that the Wilson loop insertion 
is a Lorentz-covariant expression. 

Recall that in the well-understood QED^, Polyakov j|] demonstrated that the expectation value of a Wilson loop 
operator can be written in terms of a CG with the insertion of the operator, 

(M/)^ei''(-^^), v{x,C)= I dV(2/)r^ = -/ rfV^(y)^^^'''a,^^ (21) 

JSc F ~ V\ Jsc ^ F ^ 2/1 

where the factor of 1/2 appearing in the exponential is due to the external quarks being in the fundamental repre- 
sentation; r\(x, C) is the solid angle subtended by the loop C at the point x; and Sc is an arbitrary surface bounded 
by the counter C. For simplicity, will be assumed to lie in the {xi,X2) plane. In the case of infinitely large 
Wilson loops, the solid angle is 2tt below the surface, and changes discontinuously across the loop to — 27r. Physically 
r]{x, C) is the magnetic scalar potential due to a dipole layer of unit strength on the sheet Ep- As expected, formula 
( pT| ) can be expressed exclusively in terms of the 3d Greens function which is the only relevant clement of the CG 
representation. Also note that, if the fractional monopole charges of strength 1/iV were present in the system, an 
additional factor of 1/iV would appear in the expression for ■q{x,C) in Eq.(^). From Eq. (|l|), Ur^ can be easily 
computed, 

□?7 47resc(a;i,a:2)'5'(x3), 7^ = 27resc (a^i, a:2)sgn(x3), (22) 

where QY.ci^^i^'i) unity within the surface Sc and zero outside. 

As argued above, the transition from 3d to 4d is realized by the replacement 

^ ^ 1 -> / di, (23) 



47r|x-y| 27r2(x^ - y^)2 ' 

where the numerical factors, 47r and 27r^, are the volumes of the two- and three- dimensional unit spheres respectively, 
and appear so that the Green functions arc properly normalized, OG{xy) = —5{x — y). Furthermore, an integration 
over all time-slices must be included to complete the transition to 4d. Carrying out the replacements ( p3|) in Eq. (|21 
the following expression for ri{x,C) written in the covariant 4 dimensional form is arrived at. 



-/ d^<j^,{y)(fdx,e'^''^^dx , ^ „ , 



r7(x,C) = — / d^a^^iy) i dx^e^-^'^dx- r^, (24) 
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here the integration over time-sHces J dt for a static particle was replaced by the integration along an arbitrary 
trajectory of a particle, e^^^dt e^^^'^'^ § dx„. In the static limit, ( |24| ) is reduced to expression ( pT] ) as it must 
by construction. It is interesting to note that Eq. ([2^ ) for 77 (with appropriate normalization) is formally similar 
to the expression for the linking number of a closed oriented surface d'^a^i, and closed oriented curve § dxa in four 
dimensions. However, the surface Ec, bounded by the Wilson loop, and a trajectory of the particle are not closed 
manifolds; nevertheless rj can be interpreted as measuring a solid angle, in addition ^ has an integer-number- 
discontinuity when a particle crosses the surface Ec" analogous to the Zd case (^l|) , where for a very large Wilson loop 
^ = ±1 just above and just below the surface. 

From Eq. (p^), U-q can be calculated for the Wilson loop in the [xi, X2) plane with the following result, 

□77=1(2^2)/" d'a^^i^y) i dy„e^^-^-dx5\x-y)^AT,Q^{x^,X2) I dyxe^'^dj^x-y), A,(7 = 0,3, (25) 



and is the final expression for the source of a Wilson loop in four dimensions. For the static case, the integration 
over t can be carried out and the above expression reduces to the id result found in Eq. ( p2[ ) . For very large Wilson 
surfaces, Eq. (E|) suggests the following simplified expression for 77, 



7? = Q^{xi,X2) J dyxe^^da \n{x - yf , A, fx = 0, 3. (26) 

Indeed, applying the operator □ to Eq. (p6|), and taking into account the relation □ ln(a; — yY — 47r5^(x — ?/), 
reproduces Eq. (p5|). Equation ( p5| ) explicitly shows that r\l(2-K) has a discontinuity whenever a particle crosses the 
surface. This is precisely the property of the Wilson loop insertion. Drawing attention to the behaviour in the SG 
picture, Eq. (^, the Wilson loop insertion corresponds to the shift ^ (/)-|-77 in the flavor singlet term of the potential 
(|): £:cos(^(0 + 77-6l + 27r7i) + 27rfcy 

The above discussions should have convinced the reader that although the effective SG action for low energy QCD 
does not contain the fundamental gauge degrees of freedom, it is possible to recover the Wilson loop insertion operator 
by inserting the appropriate source term. The two key points in making such a correspondence were: 

(i) Q^''^ was identified with the physical magnetic charges; consequently, the singlet combination of fields, was 
identified with the physical magnetic scalar potential 

(ii) Confinement is realized in this system through the condensation of the magnetic charges Q'^'^^ such that ^ 0. 
In this case the insertion of the Wilson loop in the monopole plasma leads to a shift (/>+7] of the corresponding 
magnetic scalar potential as demonstrated above (pO[). 

One further point is that since the species Q}^^ have fractional charge of strength 1/A^c, this additional factor should 
be introduced in front of Eq. (p^. 

Even though the insertion of the Wilson loop operator in the CG representation is now understood, the calculation 
of its VEV is much more difficult in comparison with Polyakov's model [0. The fundamental difference is related 
to the fact that in Polyakov's case the weak coupling regime is justified, and the classical action dominates over the 
contributions from one loop fluctuations. In the present case the only relevant dimensional parameter is the vacuum 
energy, E ^ ^qcd, all other parameters are expressed in terms of E and are of the same order of magnitude; as such, 
there are no small parameters in the problem. Nevertheless, the VEV of the Wilson loop can be estimated in the 
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semi-classical approximation analogous to Ref. ||l|. In this case for each given time-slice, the problem is effectively a 
3d problem described in [|l| where the area law has been demonstrated {W) ^ exp {—S) with calculable string tension. 
In our Ad case, it is also expected that (W) demonstrates the area law in agreement with our earlier calculation of 
the magnetization, {A4) ^ 0; however, an explicit calculation is still lacking. 

To conclude this section, we note that a second important source appears in this system due to the existence of 
domain wall solutions to the equations of motion derived from (^. These domain walls interpolate between vacuum 
states labeled by the parameter k, and were discussed in some detail in fl^ . Similar domain walls are known to 
exist in supersymmetric models, and are reviewed quite thoroughly in ||]. Recently, Witten conjectured |^ that in 
the large N limit the domain walls connecting two vacua labeled by k and A; -I- 1 appear to be object that are not 
solitons, from the string viewpoint, but rather look like Z?-branes on which the string should be able to end. Such 
D-brane-like domain walls have recently appeared in |^ in the context of the AdS/CFT correspondence. In fact, our 
CG picture seems to indicate that such a phenomenon takes place in QCD as well. Indeed, a distinguishable property 
of Z?-branes is that in large-A^ limit its tension is ^ N. Domain walls described in [0 have exactly this property^. 
Additional support in favor of this identification comes from our demonstration of the condensation of fractionally 
charged ^ a/N magnetic particles; which implies that the electric charges in the system could also be fractional]] 
and the domain walls support excitations that carry electric charges 1/N. Consequently, the chromo-electric flux 
contained in the open string of the corresponding (non-critical) string theory can end on the QCD D-brane. 



V. CONCLUSION 

The most important "formal result" of this paper is given by Eqs. p^ ; which we claim is the new representation 
of low energy QCD that is appropriate for the analysis of the vacuum structure. The correspondence between 
these two representations is discussed at length in the text. The most important physical (as opposite to the "formal") 
result of the paper is formulated as a conjecture described at the end of Section III in which the charges, Qf'\ from 
the statistical ensemble ( |lT| , are identified with the instanton-quarks 0. If this conjecture is proven correct, 
enormous progress in our understanding of QCD will have been made. For example, it would fill the missing element 



of the well-developed instanton picture 1 22 1 to include the electric confinement of quarks as a natural consequence of 
the same well-known BPST-instantons psf , which have been under intensive study since the 70s. Indeed, as 
discussed in the text, our particles (which are the instanton-quarks, according to the conjecture) gain the magnetic 
charges which condense, (Ai) ^ 0. Therefore, the standard 't Hooft-Mandelstam picture |^ for the confinement would 
occur due to the instantons. 

Various arguments which support this conjecture have been given in the bulk of the text; however, one further 



^ To be more precise the wall surface tension a ~ ^qcd^^^^ + it^q^qcd^ ^ S'Hd in the chiral limit becomes a ~ ^qcd^^^^ y 

however, in the large- A'^ limit a ~ A'^ as expected. 
^ This can easily be seen from Eq.(^) by making a shift 6 6 + 2tt (which does not change the physics) and realizing that 

the system supports excitations with fractional electric charges Q M + N along with fractional magnetic charges M even for 

e = Q{mod2n). 
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reason why we believe this conjecture could be correct will now be presented. Our CG representation implies that the 
number of integrations dl^x^^^ in ([lT| ) for each configuration is exactly equal to the number, A/q, of fractionally charged 
particles, (5*^°^ ~ 1/p, for that configuration. Due to the fact that the total charge of the configuration is integer, the 
number A/q must be proportional to the denominator of the charge, i.e ~ p; while, the fractional charge Q*^*^-* ^ 1/p 
in the CG representation appears due to the 6/p dependence in the effective Lagrangian approach (^. Consequently, 
as long as there is a 6*/^ dependence in the SG representation, the number Mq must be an integer multiple of p, and 
the total dimensionality of integrations over d'^x^f^ is 4Mo = Apk. This is exactly the measure in the multi-instanton 
background with p — C2(G)! The quadratic Casimir operator C2(G) appears here because in SUSY theories this 
number determines the djC^iG) dependence. 

Indeed, it is well known since the work in that the B dependence of the gluino condensate {TrXX)g in SYM 
with arbitrary gauge group G is, 

{TrXX)e,k = {TrXX)e=o exp (i ^^^) , fc = 0, 1, ^2(6) - 1 (27) 



C2(G) 

where C2(G), the quadratic Casimir operator, becomes equal to the dual Coxeter number when the longest root 
vectors is normalized to have length one. In particular, C2{SU{N)) = iV, G2(S'0(iV)) = iV-2, G2(5p(2iV)) ^ N + 1. 
Formula ( p7| ) implies that for each given 9 there exist G2(G) degenerate vacuum states (which corresponds to the 
Witten index) for which (TrXX) differs by a phase factor of exp (^i ^Jq-j ^ ■ Then the 6 evolution from = to = 27r 
according to Eq. (^7|) simply renumbers these states in a cyclic way. The easiest way to understand this result is 
(roughly speaking) to count the number of gluino zero modes in the one instanton background (it is equal to 2G2(G)) 
such that the correlation function {TrXX, TrXX, ...)e^k ~ exp(i0), including G2(G) insertions of the operator TrXX is 
not zero. Formula (27) then follows from the consideration of this correlation function (for a complete analysis see the 
original paper [^). The number of bosonic zero modes is twice as much and equals to 4c2(G) ||2l[. The generalization 
of this result for an arbitrary number, k, of instantons is also well known and is given by the formula 4G2(G)A: 
These well known results are reviewed here to place emphasis on the one-to- one correspondence between 9/C2{G) 
dependence of the physical parameters and the number of zero modes 4G2(G)fc in the k instanton background. But 
as argued above, the total dimensionality of integrations over d'^xf''^ in our formula ( [ll| ) exactly equals to this number 
4Mo = 4G2(G)fc, if the 9/C2{G) dependence in non-supersymmetric gluodynamics remains the same as in SUSY 
theories. Soft supersymmetry breaking analysis strongly suggests that this is the case g. Therefore, our conjecture 
can be reformulated in the following way: if 9/G2{G) dependence in non-supcrsymmctric gluodynamics remains the 
same as in SYM, our particles ( pi]) could be nothing but the instanton-quarks suspected long ago Q]. The recent 
paper |3l|| , where it was demonstrated that in SYM the instanton-quarks carry magnetic charges and saturate the 
gluino condensate, also supports our picture. 

As a final concluding remark, at the intuitive level there seems to be a close relationship between our CG represen- 
tation in terms of Abelian monopoles, Q^^\ and the "Abelian projection" approach |23 2^j2^ as well as the "periodic 
instanton" analysis |2^. We also feel that there is a close connection with the work of in the description of 
infra-red QCD physics. Unfortunately, an explicit realization of these correspondences is still missing at the moment; 
however, the search for such mappings will be the focus of future studies. 
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APPENDIX A: ei-DEPENDENCE IN THE MASSLESS LIMIT 

In this appendix the 0-dependence of the partition function is proven to vanish in the limit of zero quark masses. 
This is an absolutely trivial and well-known result of QCD. It is easily understood from the effective Lagrangian 
approach (|^), where the 6 dependence can be eliminated in the chiral limit by the shift (/> — > (j) — 6. However, it is 
instructive to understand this 9 independence of Z by direct analysis of the statistical ensemble ( [l2| ) with nonzero 9 
within the CG representation. It will be demonstrated that the 9 independence is recovered in the chiral limit as a 
consequence of the long range Coulomb interactions. Neglecting the interactions will result in explicit 9 dependence 
of the partition function, which is clearly unwanted. The first point to notice is that in this limit only configuration 
which contain the species Q^*^-* are allowed. This occurs since the mass dependence in the partition function appears 
as (TOa/2)*^°, for a 7^ 0, where Ma is the number of charges of species Q'^"'. Therefore in the massless limit the 
partition function reduces to that of a single charge species 

For 7^ it is convenient to label the configurations by the number of positively, ni, and number of negatively, 
n2 charged particles (note that for 9 the ni and 712 need not be equal). In this case the total charge of the 
configuration^ = X^f^i QI"'' = '^1 ~ ^2, while the total number of pseudo-particles Mq = ni +n2 is fixed and the 
corresponding contribution, z''^'°'^\ to the partition function (|l2|) is, 

^(Mo,Q)(5)) ^ exp(-iQ6i)exp(-i^(*^«'0)) 

where, without losing generality, only the first branch with fc = appears here; the superscript (0) on Q has been 
removed since only a single charge species is present; E is the fugacity of the Coulomb gas; the combinatorical factor 
^-^1° 1^ I has been introduced for the correct counting of ni identical particles with charge (-I-) and 712 identical particle 
with charge (— ); and, 

g EE Q^°) = m - 712, Mo = ni + n2, Qi = ±1, i = l,...Mo, (29) 

For the special case of vanishing 0, the neutrality condition would require rti — ?i2 in the thermo-dynamical limit. Now, 
if the Coulomb interaction could be ignored, then ( ^ ) yields a free energy F{9) ~ — logZ(6') which is a 9- dependent 
expression. Indeed, the partition function for a noninteracting gas with a fixed total charge Q, 

Z'^iS) = exp(— iQ0) exp(— F^) = Z^^'<"QH9) -> / dMo Z^^''°^'^\9) (30) 

Mo 



In this Appendix, without losing generality, the charges Qi are rescaled to be ±1. 
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can be obtained by the steepest descent method with respect to Mq. If the size of the system, L, is large, the 
noninteracting case reduces to, 

exp{-F'"') = / dMo exp {Mo\n{L^E)) x cxp <^ Mn 1 I I In 1 I ^ 

^exp(Mo) exp(--|J-l (31) 



where the formula ln(fc!) ~ kink — k has been used, and Mq is the saddle point determined from the equation, 



ln{L^E)^llnl^^^l^^, ^ « 1, Mo^2L'E (32) 

In this approximation of noninteracting particles, the partition function is evaluates to, 

Z{9) =exp{~F{e)) ^ exp{2L'^E) J dg cxp(-i6lQ) exp |-^^^| - exp |2i;i'^ ^1 - y ^ | (33) 
Consequently, the free energy is given by, 

F{e) ~ -2EL^ fi-^\ 61 < 1 (34) 



2 , 

The evaluation of Z appearing in (|3^) reveals that the essential terms are those which contain the charges, 

IQI ^ EL^e, M ~ e (35) 

Mq 

and the approximation Q/Mq ^ 1 can be justified for small 9-^1 which will be assumed henceforth. Therefore, as 
expected, the free energy F{d) for the gas of the noninteracting particles explicitly depends on 9. 

The discussion now returns to the full long-range interacting Coulomb gas ( |2^ ) . This problem is in all respects very 
similar to Polyakov's model 0|. It is known from [ p3[ that the Coulomb interactions in a plasma cannot be ignored 
in this case. Hence, formula ( ^ ) is certainly misleading. In order to derive some estimation for the interacting case, 
once again the system is placed into a box of size L. In this case, excessive charge is deposited on the walls of the 
box, and the free energy is that of a neutral Debye plasma plus the Coulomb energy, similar to what happens in the 
well- understood QED^ case |33|] , 

exp(-F'^) - / dMo exp (A/oln(L E)) x exp <^ (In 1 I lln ^) ^ Jjl) (3^) 

The partition function with a fixed total charge Q can now be estimated as in ( |3l| , |3^ ), with the only difference 
occurring in the last term which is due to the Coulomb interaction. This term plays a key role in the following 
integration over all charges J dQ, 



Z(0) = exp(-i^(0)) ~ exp(2L4i^) j dQ exp{-i9Q) exp exp|-|^| 
Taking into account the extra factor exp(— the charges of the essential configurations behave as, 



(37) 



\Q\^L^9, -^^^^OatL^oo (38) 
Ivlp L'' 



consequently. 
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F{e)^-2EL^ + L^e^ = -2EL*i^l~—j (39) 

Thus, as L — > oo, the free energy F(6) in QCDi in the chiral Umit does not depend on 0, which is precisely how it 
should behave and is exactly what is expected from the Sine-Gordon representation (|^) of the partition function. It 
should be emphasized once more that this result is a direct consequence of the strong Coulomb interactions in the 
system. 



APPENDIX B: LARGE GAUGE TRANSFORMATIONS AND MONOPOLE CHARGES 



In this Appendix the large gauge transformations generated by the operator iV, as discussed in the text (13), will 
be elaborated on. It will be demonstrated that such a transformation is equivalent to the identity operator. In 
the following analysis it will be assumed, along the lines of 't Hooft [||, that formula (^3|) has much more general 
applicability and is not necessarily constrained to the weak coupling limit of the Georgi-Glashow model where it was 
originally derived. Due to knowledge of the ^-dependence in the effective low energy potential on both sides of Eq. 
(|l3|), the magnetic properties of the relevant fields describing the large distance physics will be identified. 

The strategy is to begin with the QCD effective anomalous potential before the glue-ball degrees of freedom are 
integrated out. Next, a large gauge transformation will be performed on this effective potential. It is clear that this 
operation is equivalent to applying the identity operator; however, in the course of the following calculations, the 
magnetic and electric properties of the fields making up the effective potential will be identified.. 

Recall that the effective potential is defined as the Legendre transform of the generating functional for zero mo- 



mentum correlation functions of the marginal operators G^^G^^,, G^^G^jy and rn^'*^'*, see [|12[, 17| for details. It is 



a function of the effective zero momentum fields h, h which describe the VEV's of the composite complex fields H, H, 

J dxh = I [ dx h\ , f dxh — / f dx H 
where, 

2 V 4as q in J 2 \ q An ) ' 

The potential is also function of the unitary matrix C/y corresponding to the 75 phases of the chiral condensate: 
= — K^'l'I'j?)! Uij. As explained in Section 2 the corresponding Lagrangian reproduces the chiral and confor- 
mal anomalies and also corresponds to the following behavior of the (2^)**^ derivatives of the vacuum energy in pure 
gluodynamics, 

:i2fer^ //IN n 2k / , X 2fe 



vac (0) 



Q02k 



nrf'x.(Q(xi)...Q(x2fe))^ f^j , (41) 

where Q ^ G^,^Gfj,^. This is a consequence of the solution of the U{1) problem when Veneziano ghosts saturate all 
relevant correlation fimctions (|4^) p^ ]. 

Our starting point is then the effective QCD potential W{h, U) which in Minkowski space takes the form p^ , 

iV /. . h r . h 



hcxo g— 1 

-iVW{h,U) 



e 

n= — co k=0 



2^ Z^exp l/ilog— + /ilog 

— _rv-i t— n ^ ^ 



2eE 



I , q 6 — ilosDetU + 2Tm\ h — h , x , , , 

+ iTrV{k+- + -VTt(MU + h.c.))- , (42) 

p 2tt J 2i 2 
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where M = diag(TOi| ($'^'') |) , V is the volume of the system and the complex fields /i, h are defined as in Eq. (^) . 

The anomalous effective potential ( ^2[) contains both the light chiral fields U and heavy "glue-ball" fields h , h, and 
is therefore not an effective potential in the Wilsonian sense. On the other hand, only the light degrees of freedom, 
described by the fields U are relevant for the low energy physics. An effective potential for the U , U fields can be 
obtained by integrating out the ft., h fields in Eq.(^). In SUSY models the transition from the effective potential 
for the U, h fields to the effective potential for the U fields, by integrating out the h, h fields, is analogous to the 
transition from effective Lagrangian for SQCD to the Affleck-Dine-Seiberg low energy effective Lagrangian. 
The corresponding derivation for QCD was described earlier in ||l^,|l^. The goal here is to study the properties of 
the various fields from ) under large gauge transformations. In order to gain an understanding of what changes 
this large gauge transformation produces, the complex h, h- fields are parameterized by introducing two real fields p 
and Lu, 

h^2EeP+''^ , h:^2EeP-''^ . (43) 
The summation over the integers n in Eq. (B3) then enforces the quantization rule due to the Poisson formula, 



^expr2 

n. ^ 



j,inlv^^—^]=Y,S(-VEePsmu;-m\, (44) 



This reflects quantization of the topological charge Q = {j; VEe^ simSj which must be an integer, m, in the original 
theory. Applying a large gauge transformation, T ~ e^^'^^ , then implies that m in Eq. (Q) undergoes the shift 
m — > m — 1. Indeed, by definition 

n m 

where the \9) state is the eigenstate of the large gauge transformation operator, T, with eigenvalue e^*^. Consequently, 
when a large gauge transformation is performed the topological charge of the configuration becomes shifted Q Q + 1. 
This corresponds to the replacement m ^ rn — 1 in Eq. (^^. It is quite clear that such a replacement does not 
change any physical content of the theory due to the summation over m in Eq. (jij). Therefore, the large gauge 
transformation is equivalent to the identity operation as expected. However, the transformation properties of the 
individual fields under such an operation is what needs to be understood. Hence, proceeding further, and using (|44|), 
Eq.(^2[) (with the mass term omitted) is placed in the form, 

+ 00 q—1 



^-iVW 



TH — — oo — 

where, 

9k = 9-i\ogDetU + 2n-k . (46) 

q 

To resolve the constraint imposed by the presence of the (5-function in Eqs.(^^,(|45|), a Lagrange multiplier field, $, 
is introduced, 



S{VE - sinuj — m) oc J D ^ exp ^i^V E sin uj — 



P 



m (47) 
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Wick rotating to Euclidean space by the substitution iV V, from Eqs.(p4|),p7|) one obtains, 

^ ( +00 q-l 

W{U, p, Lu, ^) = -—\ogl ^ ^eyij>[-~VEeP {{p-l)cosu) -^sinuj) 



m= — 00 fc=0 

ilogDetU + 2'Kk- + -UJ --^ -£— \ . (48) 
q q q ) VE\ J 

The last term has been introduced to regularize the infinite sum over the integers m, and the limit e ^ will be 
taken at the end, but before taking the thermodynamic limit 1^ ^ 00. Note that Eq.(^8|) satisfies the condition 
W{io + 27r) ~ W(lo) as it should since lo is an angular variable (^). Note also that Eq. ( ^8| ) has an explicit 27r 
periodicity in d. As mentioned above, the large gauge transformations lead to the shift m m — 1; therefore, the 
extra term which appears as a result of the action of the operator T is the phase proportional to m in Eq. ( ^ ) , 

y ^ g-«(e-»logr>et(7+27rfc| + £^-£<I.)_ ^^^g^ 

However, this operation must be the identity; consequently, this phase must be unity at the very end of the calculations. 
Rather than imposing this condition immediately, the transformation can now be used to understand how this identity 
is realized in terms of the contributions from different fields and/or vacuum charges at infinity. To answer this question 
the VEV's for each field entering Eq. ( ^ ) must be calculated in the thermodynamic limit. The thermodynamic limit 
1/ — > 00 of the potential (^ ) can be obtained by making use of the Jacobi identities, 

03(u,x) = ^ y exp <^ I = y exp[-?27r2^ + 2iZH (50) 
\/t:x ^-^ I X I 

k— — oo I — — 00 

which allows Eq.(p8|) to be rewritten as, 

-00 1 



W{U,p,uj,^) = -—log I ^exp[-yi;e''{(p- l)cosw-$sinw} 



- — ( e -ilogDetU + 2TTk- + -UJ --^ -2mi 
Ae \ q q q 



(51) 



Here, an irrelevant overall infinite factor ~ e"^/^ has been ignored. Eq.(51) is the final expression for the effective 
potential W ^ which is suitable for analysis of the thermodynamic limit. Such an analysis was carried out in |p7| where 
it was shown that the coordinate-independent solutions, corresponding to the vacuum states, take the following form, 

{uj)i^ ~-{e-i\ogDetU) + 2TT-r~2nk , (/?> = , ($> = , r = 0,±l,..., k = 0,1, . . . , q - 1 
P P 

(52) 



The combination 27r^r — 27rfc which enters this equation can be represented as 



2TT-r ~ 2nk ^ — I - 2TTr' r, r' = 0, ±1, . . . , / = 0, 1, . . . ,p - 1 (53) 
P P 

Eqs.( p2t^3| ) show that there are p ^ Nc physically distinct solutions of the equation of motion for the oj field, while the 
series over the integers r' in Eq. (p3[) simply reflects the angular character of the cu variable, and is therefore irrelevant. 



Substituting Eq.(52) into Eq. ([49|) , the phase which appears due to the insertion of the large gauge transformation 
operator is seen to be identically zero (mod 2?™) as expected. What is important now is the possibility to interpret 
this result in terms of the speciflc superposition of the following operations: 
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(i) e~'°s^etc^ = g*0(2:) local field has already been identified with the magnetic scalar potential. Therefore, 
a nonzero vacuum expectation value for (e"^*-^^) implies the monopole condensation and nonzero magnetization 
as discussed in section IV. 

(ii) e~'?". The local field Lo was defined as the phase of the fields h, h ( ^3| ) which describe the vacuum expectation 
value of the composite fields H,H (|4C|). A nonzero vacuum expectation value for this operator (e"*'") also 
implies that the condensation of monopoles is related somehow to the glue-ball degrees of freedom iJ, H (|4^) ; 
however, a precise statement cannot be made because the CG representation (|l^) was derived from SG effective 
Lagrangian (0) after the integrating out all glue-ball fields h, h. 

(iii) e"'*-^^^'^'^ 9''. This operator is not related to a local field, but rather, is related to the background charge at 
infinity just like in the well-known example of the 2d Schwinger model. 

To conclude: In section III it was demonstrated that the charges Q^^^ from our system ([l^ ) can be identified with 
the magnetic charges (p5|), and the field <j) can be identified with the magnetic potential (p^). In this Appendix 
the next step in analysis of the CG-SG correspondence has been carried out. The large gauge transformation ( p^ ) 
was performed on the potential and it was demonstrated that the expected identity is realized as the superposition 
(combination) of three different non-zero contributions related to the condensation of: 1) the quark fields, 2) the gluon 
fields and 3) the background charges at infinity as explained above. 
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